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1.  INTRODUCTION 

The  relative  motion  between  a spacecraft 
and  the  plasma  surrounding  it  in  the  ionosphere 
and  magnetosphere  may  be  investigated  directly 
using  planar  ion  sensors  to  examine  the  flow  of 
plasma  past  the  spacecraft.  If  the  plasma  is 
taken  to  be  corotating  with  the  earth,  as  in  the 
case  below  about  400  km  at  middle  and  low  lati- 
tudes then  all  relative  motion  is  due  to  the 
spacecraft  velocity.  This  is  usually  known 
independently  and  the  planar  ion  sensors  may 
then  be  used  to  determine  spacecraft  attitude. 
(Sagalyn  and  Smiddy,  1969a) 

If  the  ambient  plasma  has  a bulk,  drift 
motion  In  a frame  of  reference  that  rotates  with 
the  earth  then  the  situation  is  more  complex. 

In  this  case  if  spacecraft  attitude,  velocity, 
and  plasma  density  are  known  in  the  corotating 
frame  of  reference  then  the  ion  sensor  data  may 
be  used  to  determine  the  magnitude  and  direction 
of  the  plasma  flow  in  this  frame  of  reference. 
These  conditions  hold  generally  above  about  60 
latitude  and  in  these  regions  increase  with 
height  from  about  1000  km  up. 

In  this  paper  we  outline  the  steps  required 
for  processing  data  from  an  array  of  planar  ion 
sensors  mounted  on  a spin  stabilized  satellite. 
We  shall  assume  that  ambient  plasma  density, 
spacecraft  velocity  and  attitude  in  an  external 
frame  of  reference  are  all  known  independently, 
and  proceed  to  derive  the  bulk  motion  of  the 
ambient  plasma  in  this  same  frame  of  reference. 
The  determination  of  spacecraft  attitude  in  the 
zero  plasma  motion  regime  is  a simplified  form 
of  these  results. 

2.  DETERMINATION  OF  THE  PLASMA  FLOW  DIRECTION 

In  this  section  we  derive  expressions  for 
the  Pitch  and  Yaw  of  the  plasma  flow  relative  to 
the  spacecraft  coordinate  system,  using  the 
currents  flowing  to  an  array  of  the  positive  ion 
sensors.  The  spacecraft  velocity  relative  to 
the  plasma  contains  components  due  to  the  space- 
craft velocity  and  to  any  drift  or  'bulk'  mo- 
tion of  the  plasma.  To  evaluate  the  component 
due  to  plasma  motion  it  Is  necessary  to  knew 


the  spacecraft  attitude  and  velocity  in  an 
external  coordinate  system,  and  this  is  dealt 
with  in  Section  3.  First,  however,  we  must 
derive  the  Pitch  and  Yaw  angles  of  the  plasma 
flow  in  the  spacecraft  geometric  coordinates. 
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Figure  1.  Geometry  of  Sensor  in  Ion  Flow 


In  Figure  1 we  show  x y z,  the  spacecraft 
coordinates;  ON  is  the  normal  to  the  planar 
sensor  aperture,  and  OQ  (q)  is  the  relative 
plasma  flow  vector  with  Pitch  (Et  + Pq)  and  Yaw 
Yq, Angle  E-  is  a time  dependent  measure  of  the 
rotations',  position  of  the  spacecraft  as  it  spins 
about  the  y axis.  OV  (v)  is  the  satellite  velo- 
city vector,  which  at  time  t is  Inclined  to  the 
x z plane  at  angle  ot.  This  angle  varies  with 
time  only  when  the  true  spin  axis  of  the  vehicle 
does  not  coincide  with  a spacecraft  geometric 
axis,  and  this  situation  is  dealt  with  in  Section 
4. 

The  direction  Cosines  of  0Q  are: 

COS 

CoS  YQ/^(Et  + Pq) 

Sin  Yq 

Cos  Yq  Sin  (Et  + PQ)  (1) 
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Similarly,  tor  ON 


Cos  Y Cos  Ps 
Sin  Y 

Cos  Yg  Sin  Ps  (2) 

where  P and  Yg  are  the  Pitch  and  Yaw  of  the 
sensor  normal. 

For  unambiguous  measurement  of  Pq  and  Yq 
two  orthogonal  pairs  of  such  sensors  are  required 
and  the  actual  configuration  used  is  shown  in 
Fig.  2.  z 

Lj-^Y 


Figure  2.  Front  View  of  Ion  Sensor  Array. 


I * t Ne  Aq  Cos  0 (3) 

where : 

§ = transmission  of  the  aperture 
N = ambient  ion  density 
e * electron  charge 

q = relative  velocity  between  plasma  and  space- 
craft 

A * aperture  area 

0 = angle  between  q and  aperture  normal 
(NOQ  in  Fig.  1) 

Note  that  q may  or  may  not  contain  component 
due  to  plasma  flow. 

Using  the  direction  Cosines  of  OQ  and  ON 
from  (1)  and  (2)  to  derive  CoS  0,  we  obtain: 

r CO£ 

I - l Ne  Aq  [Cos  Yq  Cos  YSa(Pq  - Ps  + it)  + 


+ Sin  Yq  Sin  Yg  j 

(4) 

Let  Tj  = § Ne  Aq 

(5) 

Y = *Q 

(6) 

X = Pq+a-Y+?'t 

(7) 

The  lower  three  sensors,  #1,  #2,  and  #3  measure 
Yaw,  and  the  vertical  pair,  #2  and  #4  measure 
Pitch  as  the  array  faces  the  plasma  flow  vector 
once  during  each  spacecraft  spin. 

Realistically  we  must  take  account  of  any 
minor  offsets  in  the  mounting  of  such  a sensor 
array  relative  to  the  spacecraft  geometric 
axis  x y z.  If  the  centerline  of  the  array 
(line  bisecting  normals  of  sensors  #2  and  #4), 
nominally  on  the  x axis,  is  rotated  y degrees 
about  the  y axis  and  e degrees  about  the  z 
axis,  then  we  have  fully  described  the  alignment 
of  the  sensor  array  relative  to  the  spacecraft 
axes.  The  angle  between  the  normals  of  any  two 
adjacent  sensors  is  2a  so  that  the  values  of 
Pg  and  Y§  for  the  four  sensors  are  as  follows: 

TABLE  I 

Sensor  No.  1 2 3 4 


From  Equation  4 and  the  angles  given  in  Table  I 
we  obtain  the  following  relations  where; 


Rij 

- (It  - Ij)  / (Xt  + Ij) 

(8) 

r21 

= tan  a Cos  X (tan  a + e)  + tan  Y 
Cos  X - tan  Y tan  ( a + e) 

(9) 

R13 

= -tan  2a  Cos  X tan  e + tan  Y 
Cos  X - tan  Y tan  e 

(10) 

R23 

= ten  a Cos  X tan  (a  - e)  - tan  Y 
Cos  X + tan  Y tan  (a  - e) 

(11) 

h + 

1^  = 2 I2  Cos  2 a 

(12) 

Now 

let  tan  $ = tan  Y/Cos  X 

(13) 

then 

R21  = tan  a tan  (a  + e + $) 

(14) 

R23  = tan  a tan  (a  - c - $) 

(15) 

R^j  * -tan  2 a tan  (e  + i) 

(16) 

-(ar-v) 

-(2ar+c) 


-(o-Y) 


-(or-Y)  (cy+y) 

-(2a“t)  -e 


Considering  again  the  Ion  flow  to  a single 

sensor,  we  assume  that  the  mean  thermal  an(j  „ _ Sin  a Sin  (a  - XI 

(Isotropic)  velocity  of  the  Ions  is  much  less  ^4  Cos  a Cos  (a  - X)  - tan  Y tan  e (17) 

than  the  spacecraft  velocity  (Sagalyn  and  Smiddy, 

1969b)  so  that  I,  the  sensor  current  is  given  by: 
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Define  X such  that: 

R^j  + tan  2 cr  tan  e 
* R13  - t an  2 o e 

CPT 


R24  -1--X-  ■■ 

mi, 


= tan  a tan  (a  - X) 


Explicitly : 

999 

f * q + v • 2q  v (Cos  of  Cos  Yq  Cos  Pq  + 
+ Sin  ot  Sin  Yq) 

The  direction  cosines  of  the  true  plasma  flow 
are  given  by: 


In  equation  19,  X is  the  only  unknown, 


tan  (P0  + E,  - V)  = -R94 1 - 

tan  a 1 + X R, 


We  can  now  insert  value  of  X just  evaluated, 
into  equations  14  or  15,  and  16  to  determine  Yq. 
The  choice  between  R2 ^ and  3 depending  on 
which  side  of  the  centerline  the  plasma  flow 
is  located. 

3.  DERIVATION  OF  PLASMA  VELOCITY  VECTOR 


cos  a - [v  cos  ot  cos  - 

- q cos  Yq  cos  <Et  + Pq)  ] /f  (23) 

coo  p = [V  sin  at  - q sin  Yq  ] /f  (24) 

cos  M = [ V cos  Sin  Et  - 

- 9 cos  ^ sin  (Kj.  + Pq)  ] /f  (25) 

The  pitch  Pp  and  Yaw  Yp  angles  of  the  true  plasma 
flow  are  related  to  these  direction  cosines  by: 

cos  Yp  cos  (Et  + Pp)  - cos  a (26) 

sin  Pp  = cos  0 (27) 

cos  Yp  sin  (Et  + Pp)  = cos  M (28) 

Thus,  the  pitch  and  yaw  angles  can  be  determined 
in  terms  of  the  magnitudes  and  directions  of 
relative  plasma  flow  and  vehicle  velocity; 

tan  (Kt  + P ) - 

V cos  <rt  sin  E£  - q cos  Yq  sin  (Ej  + Pq)  (29) 

V cos  a t cos  Efc  - q cos  Yq  cos  (E,.  + Pq) 

tan  Yp  = 

V sin  " q sin  Yq 

£v^  cos^  <Tt  + q^  cos^  Yq  - 2Vq  cos  cos  Pq  cos  Yq j ‘ 

(30) 


Figure  3.  Definition  of  Ion  Flow  Vector. 


In  Figure  3 we  show  the  vector  triangle 
containing  the  true  plasma  flow  vector  f 
together  with  v and  q from  Section  2.  The 
angle  9,  between  v and  q is  determined  knowing 
Pq  and  Yq  from  Section  2.  The  value  of  v,  the 
spacecraft  velocity  vector  is  supplied  inde- 
pendently, and  q,  the  relative  velocity  between 
the  plasma  and  the  spacecraft  is  derived  from 
Equation  1 where  q is  now  the  only  unknown. 

In  Figure  3 we  now  have: 


2 2 

q + v - 2q  v CoS  0 


4.  THE  EFFECT  OF  SPIN  AXIS  OFFSET 

The  previous  results  have  been  derived  in 
terms  of  Instantaneous  values  of  all  variables. 

In  this  section  we  consider  the  effect  of  a space- 
craft spin  axis  displaced  from  the  geometric 
x y z axes  in  Fig.  I.  This  is  a realistic 
condition  since  in  flight  the  exact  location  of 
the  spin  axis  will  be  set  by  the  deployment  of 
antennas  and  other  items.  This  Introduces  a 
precession  of  the  spacecraft  coordinate  system 
x y z about  the  true  spin  axis , and  the  true  spin 
axis  sweeps  out  a circular  precession  conic  about 
the  spacecraft  y axis,  with  half  angle  A at  the 
spin  frequency  u . 


Let  us  now  define  a new  spacecraft  coordi- 
nate system  x'y'z'  in  which  y'  is  at  angle  A to 
the  y axis  and  for  simplicity  of  illustration 
x'  is  taken  to  coincide  with  the  x axis.  If  x' 
and  x do  not  coincide  then  a rotation  about  the 
y axis  is  required,  affecting  transformation  31. 
If  all  spacecraft  attitude  data  is  supplied  in 
the  'prime'  system  in  order  to  remove  the 
modulation  caused  by  spin  axis  offset  then 
we  must  establish  a means  of  using  these  data 
with  the  instrument  outputs  in  section  3 in 
order  to  derive  the  plasma  magnitude  and 
velocity  in  any  coordinate  system  exte*nal  to 
the  spacecraft.  The  two  spacecraft  coordinate 

systems  at  time  to  when  the  spin  axis 
lies  in  the  yz  plane,  are  related  by: 


Referring  to  Fig.  1,  the  pitch  Kt  and  yaw  at 
angles  of  the  velocity  v are  obtained  In  terms 
of  it  * (t0) , o'(t0),  0,  and  w as: 


Kj,  “ tan 


t 

) COS  K'  (t0)  + ut 


cos  b’1  (t„)  cos 8 sin  u 


+ sin  # (sin  K*  (t  ) cos  ut  - tan 


t = Tan*1  | j^Sin2  (fc'(to)  + ut)  + Tan2  o'(t0) 
j^Cos  K'(t0)  Sin  # Sin  ut  + Tan  <r'(t0)  Cos  8 + 
+ Sin  K'(t0)  Sin  8 Cos  utl/ 


0 Cos  A -Sin  A 


0 Sin  A Cos  A 


where  the  transformation  matrix  is  orthogonal. 

To  transform  velocity  v*  as  supplied  by 
the  attitude  system  under  these  operating 
conditions,  to  v in  the  x y z system: 


5.  APPENDIX 

The  transformation  properties  of  four 
useful  coordinate  systems  are  listed  in  Table  II. 
The  systems  are: 

Earth  Centered  Inertial #....e 

Orbital  Coordinates 

Local  Vertical. .....1 


Spacecraft  Coordinates. 


COORDINATE  TRANSFORMS 


Cos  wt  0 -Sin  ut  Vxto 


0 Cos  A Sin  A 


0 -Sin  A Cos  A I Sin  wt  0 Cos  wt 


Explicity  the  velocity  vector  V(t)  at  time  t in 
spacecraft  x y z coordinates  can  be  obtained 
from  V'to  by: 

vx  (t) 

Vy  (t)  = 

Vz  (t) 

Vx(t0)  coswt  - V£(t0)  sirovt 

v;<t0)  sin#  sinut  + Vy(tc)  cos#  + Vz(t„)  sin#  cos  t 
Vx(t0)  cos#  sinut  - Vy(t0)  sin#  + Vz(tQ)  cos#  cos  t 
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f 


*:'*v 


For  example: 

e - ETS,  o - De 


where  lower  case  Indicates  3x1  column  matrices,  and  capital  denotes  3x3  square  matrices. 


The  transforms  are 

as  follow: 

A * 

0 

cos  » cos 

cos  »g  sin  t Q cos  * 

■ sin  i sin 

o ° 

COS  v 

0 sin  sin  <PQ  + sin  o0  cos 

-sin  i 

0 

cos  i 0 COS  <PQ 

cos 

o sln  o 

-sin  v cos  <!> n 
o o 

-sin  *0  sin  cos  *0 

- cos  c0  sin  •fig 

-sin 

sin  sin  <PQ  + cos  »0  cos  v>o 

Aj  * Aq  with  subscripts  o replaced  by  1. 

B - cos  8 o 

o 1 

-sin®  o 


cos  8 sin  1 
cos  i 

-sin  8 sin  1 


D - 

- COS 

i sin  12 

cos  i cos  12 

sin  i 

- sin 

i sin  12 

sin  i cos  12 

cos  i 

_ - cos 

(2 

- sin  12 

o 

- cos  8 cos  1 sln!2  - sin  8 cos  12 

- sin  1 sin  S3 

sin  9 cos  1 sm  12  - cos  6 cos  SI 


cos  8 col  1 cos  12  - sin  9 sin  12 
sin  1 cos  12 

-sin#  cos  1 cos  12  - cos  6 sin  12 


sin  6' 
o 

cos  8 


E - AQD  - AjC 

where  superscript  T denotes  transpose  matrix, 

8 » latitude 
12  » right  ascension 
•C  - roll 
^ “ yaw 
* “ pitch 
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13.  ABSTRACT 


— > I.ow-energy  plasma  motion  in  the  ionospheric -magnetispheric  region 
mav  be  measured  directly  using  an  array  of  planar  ion  sensors  on  a spin- 
stabilized  satellite.  In  the  equatorial,  low-altitude  (below  400  km).,Jlzero 
plasma  motion*1  region,  the  relative  plasma  flow  is  due  to  spacecraft  velocity 
onlv;  the  sensors  give  information  concerning  the  attitude  of  the  spacecraft  in  the 
plasma.  In  the  polar  regions  or  at  higher  altitudes,  plasma  flow  is  important, 
and  the  attitude  and  ephemeris  data  of  the  spacecraft  in  an  external  system  are 
necessary  for  the  determination  of  plasma  flow  in  this  system.  We  outline  a 
sensor  system  that  measures  the  How  of  lovv-energv  plasma  and  proceed  to 
define  all  the  additional  information  required  to  operate  in  the  moving  plasma 
region,  together  with  the  coordinate  transforms  this  involves.  Corrections  are 
derived  for  three  factors  that  arise'  in  the  practical  application  of  such  an 
instrument:  1)  the  offset  between  spacecraft  geometric  axes  and  the  actual 
alignment  of  the  sensors  as  mounted;  2)  the  offset  between  spacecraft  geometric 
axes  and  the  spin -axis  location  achieved  in  flight.;  and  3)  the  separation  of  the 
spacecraft  velocity  vector  from  the  spacecraft  spin  plane. 


KKYWOltIJS:  Plasma  bulk  motion,  Ion  drift,  Polar  wind,  Ion  sensors. 
Spacecraft  altitude,  Spacecraft  coordinate  t ransforms 
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